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Abstract 



A number of new definite integrals involving Bessel functions are pre- 
sented. These have been derived by finding new integral representations for 
the product of two Bessel functions of different order and argument in terms 
of the generahzed hypergeometric function with subsequent reduction to spe- 
cial cases. Connection is made with Weber's second exponential integral and 
Laplace transforms of products of three Bessel functions. 
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1. Introduction 

The aim of this work is to derive a number of infinite integrals involving 
Bessel functions which appear to be new. The approach is, beginning with 
an expression for the product of two Bessel functions as a sum of Gauss 
functions, to integrate and perform a resummation to obtain other hyper- 
geometric functions, and then to reduce these to more familiar form. In 
particular we exploit the relation between Bessel functions and the function 
oFafl]. Finally, we generahze Weber's second exponential integral by express- 
ing the Laplace transform of a product of three Bessel functions as an infinite 
series of products of modified Bessel functions. We have not aimed at com- 
plete rigor or generality; operations such as interchange of limits and Hankel 
inversion are carried out formally, i.e. without verifying that the conditions 
stated ensure their validity. In the majority of cases approriate conditions 
can be supplied by appeal to convergence and analytic continuation. 

2. On the product J^{ax)Jy{hx) 

Our starting point is the familiar expansion [2] 



V{^l + l)V{v + l)J^{ax)J,,{hx) 



(-l)"(|ax) 



2m 




m=0 



m!(/i + 1) 



m 
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By making use of the standard transformation [3], we have 
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Fi(-m, - m; i/ + 1; — ) 



a 



- (1 - ^)/^+-+2m+l ^p^^^ + ^ + 1^ ^ + ^ + ^ + 1.^+1-^) (2.2) 

(X CL 

Next, we recall that [4] 

2Fi(i/ + m + l,i/ + // + m + l;i/ + l;— ) = 



4(a/6)^ 



r(/x + 1/ + + + 1/ + 1), 

6 



/ t''+''+^"'+^KJ2t)lJ2-t)dt (2.3) 
io a 



b 

m = Q,l,2,...,Reu > -1, Re(ii + u)>-l, |i?e(-)| < 1. 

a 



By inserting (2.2) and (2.3) into (2.1), we get 

1 1 6^ 

-axn-bxni - - 

2 2 a 



r{u + l)r(/x + l)T{ix + l)J^{ax)Mbx) = 4(-ax)'^(-6x)^(l - ^r^''^'- 



(7 /""^ b 
jy ti^+-+\Fs{fi+l,i^+l,fi + iy + l;-zH')K^{2t)I,{2-t)dt (2.4) 
b Jo a 

Reu>-1, Re {/jl + u) > -1, \Re{-)\ < 1, 

a 

where z = \ax{l — ^). 

It is interesting to observe that eq. (2.4) enables us to derive the Mellin 
transform of J^{ax)J,y{bx), i.e. the Weber-Schafheitlin integral [5], in a simple 
way. Indeed from [6], 

PCX) 

- - VP-'-— ^ r(H±ii^)r(/. + i)r(. + i)r(/. + . + 1) 

i?e (// + — s) > —1, 

and [4] 
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Jo ' ^M(2^)/.(2-t)dt - ^f(^r( )r( ). 

+ ^ + - + \ + ^ + S . + 1; ^) (2.6) 

i?e(i/±// + s) > -1, |i?e(-)|<l 
we immediately obtain 



^^"^^ ^^^^ 2 ' 2 + 

(2.7) 

i?e(// + i/-s) > -1, 0<6<a 
With b ^ ib {b > 0),and x positive, eq. (2.4) becomes 

r(// + i)r{u + i)r(/x + u+ i)j^{ax)h{hx) = 

{^-axn\bxr{i+^^,r'-''^\lr- 

/ oi^3(/^+l, -— aV(l+^) V)X^(i) J,(-i)cii (2.8) 

JQ lb a 

Rev>-1, Re{n + i')>-l 
or, by writing b=ay and x=a/4 

r(/. + i)r(^ + i)r(/. + ^ + i)j,{^a')U^a'y) = (^)''+''(i + y'r^''^'- 

• / t''+-^+\F,{t, + l,u+l,t, + u+l;- — {l + yYf)K^{t)Myt)dt (2.9) 
>/o zoo 

it!e 1/ > -1, Redj^ + v) > -1. 

Also, with ^ 



1+2/^ 



r(/. + i)r(^ + i)r(/. + v + ^VA^YA^) = (i + y^)a^'^''- 

roo 

■ t^'-^''+\F3{ii + l,u + l,ii + u + l;-aH^)Kf,{t)J,{yt)dt (2.10) 

t/ 

Re u > —1, i?e(/x + u) > —l,a,y > 0. 

For i?ez/ < —1/2 we may use the Hankel inversion formula [7], and eq. (2.10) 
gives 

(ai)'^+%F3(// + l,i/ + l,/x + i/ + l;-a¥)X^(i) = r(// + l)r(i/ + l)r(// + i/ + l)- 

■ r TffMty)JAYf^,)'.(^,)''y (2.11) 

This is an interesting addition to the class of Sonine-Gegenbauer integrals 
[8]. In particular, by dividing by a^''"'' and taking the limit a — > 0, which 
may be taken under the integral sign, we have 

{-r^K,{t) ^ r(^ + . + 1) I j^r^^^^^Mty)dy (2.12) 

The four particular cases 
of eq.(2.9) are of some interest. By using 
and 

— QF3(-,-,l:- — ) = bet(x) (2.14) 
4 " ■^4' 2' 256^ ^ ^ ^ ' 

we obtain the four integral representations 

|(1 + ^2)-i/2 j^(^) cosh(^) = 1^°° Ko{t)ber{a^{l + y^)t) cos{yt)dt 

(2.15) 
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|(1 + y^)-VVo(^) sinh(^) = K,{t)hei{ayJ {I + y^)t) sm{yt)dt (2.16) 
(1 + yY^/^I,{^) cos(^) = e-%er{a^{l + y^)t)Jo{yt)dt (2.17) 
(l + y2)-i/2/o(^)sin(^) = J\-'bez{a./{U^t)Myt)dt (2.18). 



Therefore, 



Ko{t)ber{aVl) ^ I (1 + y^)-VVo(-^) cosh(-^) cos(iy)dy 

^ ^ (2.19) 

Ko{t)bei(aVl) = / (l+y^)-VVo(--— ^)sinh(-— ^)sm(iy)dy (2.20) 

V6er(av^) = |/(1 + y^)-V2j^( ^ eos(l^) Jo(ty)dy (2.21) 

V*6d(a^/t)= r|/(l + y^)-V%(i^)sin(i-^)Jo(ty)dy (2.22) 
I JO 4 1 + ?/^ 4 1 + ?/^ 

A curious formula arising from (2.4) may be mentioned. We set = |, 
replace Ii,{2H) by Poisson's integral [9] 

i?e 1/ > -- 
2 

interchange the order of integration, and then use the Laplace transform 

r e-^H-'^^' ^ + 1, ^ + ^; -ae)dt = £(^f^/5-(2^+i) sin(4^) 

Jo II 4-1/0; p 

(2.24) 

i?e i/ > -1, i?e /? > 

which is readily established by expressing the 0-^3 as its power series and 
integrating term-by-term. This leads to 

sin(aa;) J,y(6a;) = 



6 



Rev > --, \a\ > 161 > 0. 

2' ' ' ' ' 



Therefore, with x — ^ and h — ay 



T- -vryr 1 - y^)''+2 / \ sin —)du 2.26 



i?ei/>-^, < 1. 

By taking advantage of other known reductions for the o-^s functions in 
(2.4) many additional new integrals can be derived. For example from [23] 
we obtain 

/ e~^/iy(xsin 6')[sin(— — )6er2i/(2 cos 9^/ux)+cos(——)bei2u('2cos 9\/ux)]dx 
Jo 2 2 

= sec 9 sin u Jv{usni 9), 

and in particular 

poo 

/ e~^l2n(xsin ^)6ei4„(2cos 9\/ux)dx — (—l)'^sec9 sin u J2n(usin 9) 
Jo 

roo 

I e~^/2n+i(3:^sin ^)6er4n+2(2 cos 9\/ux)dx — i—lY'^^sec 9 sin u J2n+i{usm 9). 
Jo 

We conclude this section by deriving a further integral representation for the 
product J^{ax)Ji,{bx) in terms of a 0-P3, but different from (2.4). 

Let us consider eq.(2.1) with /i — 1/. According to the quadratic transfor- 
mation [10] 

= (1 + gr-^E (2.27) 
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and (2.1) becomes 



Ji,[ax)Jy[bx) = {-abx ) } ^ ■ 



m=0 

ab \2r 



r(i/ + m+ 1) 



V \a^+itl . (2.28) 

^ (m-2r)!r!r(i/ + r + 1) ^ ' 



By using [11] 



it follows that 



CC L 2 ■ 



E r)^ c(m + 2r, r) (2.29) 

m=0 r=0 m,r=0 



Ji,{ax)Jy{hx) 



^ 1 g&a: ^ (-l)"'(|xv/^^TF)'^+^''+^"^ 

, 1 abx 



We note, in passing, that (2.30) provides a quick derivation of Weber's second 
integral [12]. Indeed, by using 

^oo (n'2 -I- ^2W/2+r 2 , ^2 

I e-^-^x-+^-+V.^..(xx^^)dx ^ ^"^^+ ^V^^^^ e-"^ (2.31) 



we have 



/ xe JJax)JJhx)dx = —e ^p'' Y] 
Jo 2p ~L 



2p r!r(i/ + r + 1) 



1 ^_£it^ g (a6/4p)^+2^ 

r-=0 

1 g^+b^ ^ ,oh. , , 

= _e T^/^ _ 2.32 
We also observe that (2.30) can be obtained from the formula [13] 

2Fi(a, 6; c; x)2Fi(a, 6; c; 7/) = 

= E (")-W-(^- "M^-^)- (^^)^ (a + r, 6 + r; c + 2r; X + y - xy) (2.33) 

r=0 ^"•(,Cjr(,Cj2r 
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by applying the confluence principle [14] twice. We first replace x by x/b, y 
by y/b, and let b— > oo. This gives 

,F,{a; c; x),F,{a; c-y) = f^ (aMc - a) + ^. ^ + 2r; x + ?/) 

t^Q r\{c)r{c)2r 

(2.34). 

Next, we replace xbyx/a, ybyy/a and let a— > oo; then 

oF,{c; x) oFi(c; x)^J2 , A / ^ oi^i(c + 2r; x + y) (2.35) 

which is equivalent to the desired result. If on the right hand side of (2.30) 
we write [9] 



((l/2),TVa' + 62)-+2r 



we get (i?e u > 



[\l - t^y-h+^^ cos(xtV^iFT¥)dt, 
Jo 



(2.36) 



J,(ax)J,(bx) = —J-—(ahxy [\l-tY-'^cosixtV^t?+¥)- 

TTL [ZV + 1) Jo 

■ oFsii^ + 1' i + i' i + - '" f)dt (2.37) 



and with 6 = 1, x = a^^s^ and a = i(Vii^ + 2 + y/u^ - 2)^ 

2 2 



2 



Jo 



7rr(2i/ + l) 

z/ 1 z/ 31 
2 ^4' 2 ^4' 64 
Finally, with i/ = 1/2 (2.37) yields [23] 



oF^{v + 1, ^ + ^ 1 + r; _L(i _ t'f)dt (2.38) 



= (2/m) sm{u^a/2b) sin{u\Jh/2a). 
Further, more complex, evaluations are possible by the same procedure. 

3. The integral e-''^ Jo{l3i^)Jo{l32V^)Jo{l3zV^)dx. 

The integral in this section heading, which we denote /32, /^s) is an 
extension of the i/ = case of Weber's second integral, eq. (2.32), to which 
it reduces when one of the three parameters vanishes. 

By writing (see eq.(2.1)) 

Jo(/32v^) Jo(/33V^) = E 2Fi(-n, -n; 1; §) (3.1) 

n=0 

and [15] 

re-'^^x-JMV^)dx = ^e-^?/^"L4^), (3.2) 
Ln{x) being a Laguerre polynomial, we first have 

/(ft, A. ft) = -e-^ E iz_L(^). ,f,(_„,_„; :;|)i„(^) (3.3) 
Now [16] 

2Fi(-n, -n- 1; |) = (1 " f)'^ 2Fi(-n, n + 1; 1; -^^) 

= (l-im(§^) = -^ r(/32' + /3|-2/?2/?3COS ^)'^d^. (3.4) 
Therefore, in terms of a well known generating function [17], (3.3) becomes 

/(A, /32, /^s) = ^ / E ^^^(/52 + - 2/^2/33 COS eru^)de 

na Jo nl Aa 

= ^ / e-^^/^«/o(§^)d^ (3.5) 

where i? = (/3|+/3| — 2/32/93 cos ^)^/^. Next, by using Graf's addition theorem 
[18] 

^o(§f ) = E (2 - U4(§^)/n(§^) cos(n^) (3.6) 
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and the familiar integral representation 

Te^- ^cos(n^)d^ = nU^), (3.7) 
Jo 2a 

we finally have 

/•CO 

/ e-"^Jo{PiV^)Jo{P2^/^)MP3V^)dx 
Jo 

-1 oo 

= ^e-^JPl+(^HPl) ^(2 - Son)FM, P2, Ps) (3.8) 

where 



« n=0 



F4(3u (3., /?3) = ^"(§^)^4§f^)/n(^) (3.9) 
We conclude by sketching the evaluation of the more general integral 

POD 

/m(/3l,/52,/53) = / e-''^M/3l^/^)Jm{/32^/^)M/33^/^)dx (3.10) 

Jo 

where m is a positive integer. In the first place, (3.4) takes the form 

E ^(i)" ^F.(-n,-n -m;r„+l; ("D 

By observing that [19] 

Finn m-m^V^h-il c^+m fi + f^^ 

= M ^2n / sin^- ^ (/?2^ + (31 - 2(32^3 COS ^)"cie (3.12) 

[Zmjlnp^ Jo 

eq. (3.11) becomes 

UPM = /'sin-* t tllpl. 

ira [2171)1 a Jo nl 
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{131 + 131 - 2P2Ps cos eTL^+n{^)de. (3.13) 

Now, from [20], 

77 ^ d 

and [17] 

oo j.n 

^ 7;;7t;;tt^" = (a;i)-/2e*j^(2v^) (3.15) 

n=0 + "J- 

it follows that (where R has the same meaning as before) 

£;^(-f;)"W,(g)^e<«-./-2".(A)„. 

{R.^)-"'UR'j^%.eii>„- (3.16) 

By inserting (3.16) into (3.13) and using [21] 

{rJI)-^I^{rJI) = {^)-^{m - 1)! f;(m + n)C;"(cos 9)- 

V CX V (X ZiGL ^ Q 

Im+n{P2\f^)Im+n{P3\f^) (3-18) 

V a \ a 

and [22] 

sin'"^ ee^""' ^Cr*(cos e)d^ 

JO 

_ 7r2^-(2m + n-l)! 2a /^^/^s. 
- n!(m-l)! ^^"^^^ 
we finally obtain 

/■oo 

/ e~'^''Jo{Pi^/x)Jm{P2Vx)Jm{P3Vx)dx 

Jo 

^ (n + m)(2m + n-l)! ^ /x /?2/?3^| 



12 



In particular, by letting /^a — > 0, after dividing by (3^, 
Jo 

^ ie-/^i/^"(;^)-e-(^)™/^/„(/3,w^)L=,./,, (3.20) 
a ax a \ a ^' 

a^2a> to \nrP2^ 2a 

which is of interest in connection with formulas (39)- (42) on page 186 of 
reference [15]. 

In conclusion, we point out that the derivation of (2.24) can be extended 
to give 

roo 1 

e-^^ oi^3(/^, + 2^ -«'^')^^ = (2a)^-^r(/x)r(2i/)/5'^-2--^ J^_i(4a//5) 

(3.21) 

Therefore, the results of section 2, for example, are capable of extension 
in a variety of directions. We leave this for the future and merely quote 
one example, interesting because it contains each of the four types of Bessel 
functions: 

/■oo 

/ xJAax)h{ax)Yoix)Koix) dx = -(27ra^)-^ ln(l - a^) (3.22) 
Jo 

where < a < 1. 
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